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A new class of accelerating, exact, explicit and simple solutions of relativistic hydrodynamics is 
presented. Since these new solutions yield a finite rapidity distribution, they lead to an advanced 
estimate of the initial energy density and life-time of high energy heavy ion collisions. Accelerating 
solutions are also given for spherical expansions in arbitrary number of spatial dimensions. 
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I. INTRODUCTION 

Relativistic hydrodynamics can be successfully applied 
both at the largest and the smallest scales of physics. 
However, its equations are highly nonlinear, and only 
few exact relativistic solutions are presently known. We 
present here a recently found family of exact solutions [1] . 
These generalize the Landau-Khalatnikov solution [2-4] , 
as they are also accelerating and yield a finite rapidity 
distribution. They are simple and explicit, similarly to 
the Hwa-Bjorken solution [5, 6]. We show new 1+ 3 
dimensional solutions [1], that are reduced to the exact 
hydro solutions of refs. [7-11] in the accelerationless limit. 



II. NEW SIMPLE SOLUTIONS OF 
RELATIVISTIC HYDRODYNAMICS 

Let us denote the metric tensor by g^'^ — 
diag{l,— I,— 1,-1), the four-velocity by u^^ = 7(1, wn), 
the pressure by p, the energy density by e, the tempera- 
ture by T, the charged particle density by n and the en- 
tropy density by a. All fields depend on the coordinates 

. We denote by r the spatial coordinate in 1+1 di- 
mensions, and the radial coordinate in 1-f- d dimensions 
(where r > 0). Relativistic hydrodynamics expresses lo- 
cal momentum, energy and charge conservation: 



{g^P^u^uP)dpP, (1) 
0, (2) 
0. (3) 



From these relations entropy conservation also follows: 
du{(Ju^) = 0. Our choice of Equations of State (EoS) 
is that of an ideal gas: p — nT , e — np. For an ultra- 
relativistic ideal gas in d spatial dimensions, k = d = 
l/Cj. In the considered case, a ocT'^. 
We use the t and 77 Rindler coordinates: 

|r| < |t| : r = rsinliT^, t = ircoshr;, (4) 
|r| > |t| : r — ztrcoshrj, t = Tsinhrj. (5) 



In order to generalize the Hwa-Bjorken solution, we in- 
troduced the A > parameter, and found the following 
class of solutions: 
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The initial pressure, number density and temperature are 
denoted by pq, uq and Tq, respectively. The above forms 
are hydrodynamical solutions for the special values of A 
and K as given below. The A = 1 , |r| < |t| , iy{s) — 1 
case is the d-dimensional generalization of the well known 
inertial Hwa-Bjorken flow [5, 6], where k and d > 1 are 
arbitrary. Here ^(s) is an arbitrary positive scaling func- 
tion, given by the initial conditions, and s is a scaling 
variable, that has a vanishing co-moving derivative: 

|r|<|i| : s(t,??)= (^)'"'sinh[(A-l)r7], (10) 



|r|>|i| : s{t,v)=[^) cosh[(A-l),;], (11) 

valid for A ^ 1. If A = 1, we have s — sorj for |r| < \t\, 
while s — t/tq for |r| > \t\. 

The constants d, A and k are constrained: In 1 + 1 
dimensions, A is arbitrary but k = d = 1, except if A = 1 
when K is arbitrary. If d is arbitrary, A = 2 and k — d 
is a valid solution both for |r| < \t\ and |r| > \t\, and for 
A = 1, we have solution inside the lightcone for arbitrary 
K. Recently, T. S. Biro pointed out [12], that the A = 1/2, 
K = 1 and the A = 3/2, k = 11/3 cases are solutions for 
d — 3 dimensional expansions, which can be generalized 
easily for any d e M, choosing k = {id — l)/3. The 
function 0(A) = for A = 0, 1 or 2, while it is unity, 
(/)(A) = 1 for A — 1/2 or A = 3/2. It is introduced to 
indicate, that the solutions are dependent on space-time 



A-l 



2 



Brazilian Journal of Physics, 2006, in press 



dn,dn 



1.2- dy dy 




X=2, m=140MeV 




FIG. 1: Fluid trajectories of the A = 2 exact solution. 



rapidity 77 in the X — 1/2 and 3/2 cases only. The A 7^ 1 
solutions describe a flow with relativistic acceleration. If 
A = 2, the trajectories are uniformly accelerating, see 
Fig. 1 and ref. [1]. 



III. RAPIDITY DISTRIBUTIONS AND 
ENERGY DENSITY ESTIMATION 

For d = 1, the rapidity distribution, ^ was calculated 
in a Boltzmann approximation in ref. [1] . The freeze-out 
temperature is T{ri = 0,t = tj) = Tf. We assumed [1], 
that the freeze-out hypersurface is pseudo-orthogonal to 
m''. With a saddle-point integration, for A > 0.5, and 



m/Tf ^ 1 (where m is the particle mass), we got 
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The "width" of this distribution is 



(12) 
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The rapidity distribution has a minimum at y = 0, if 
Ay2 < (i.e. 1/2 + Tf/{Am) < A < 1), it is flat if A = 1 
or A = 1/2 + Tj/(4m), otherwise it is approximately 
Gaussian. The typical cases are plotted in Fig. 2. 

As an application, we estimate the energy density 
reached in heavy ion reactions. Let us focus on the thin 
transverse slab at mid-rapidity, just after thermalization 
(r — To), illustrated by Fig. 2 of ref. [6]. The radius 
R of this slab is estimated by the radius of the colliding 
hadrons or nuclei, and the volume is dV = {B?"iT)TQdi^Q. 
The energy content is dE — {mf)dn, where (rrit) is the 
average transverse mass at y = 0, so similarly to Bjorken, 
the initial energy density is 



[mt) 



dn 



{R?it)to drjo 



(14) 



For accelerationless, boost-invariant Hwa-Bjorken flows 
Vo = 2/: however, for our accelerating solution 



FIG. 2: Normalized rapidity distributions from the new so- 
lutions for various A, Tf and m values. Thick lines show the 
result of numerical integration, thin lines the analytic approx- 
imation from eq. (12). For A > 1 and not too big Tf it can 
be used with about 10 % error. 



we have to apply a correction factor of = 

{tj /tq)^^^ (2A — 1). Thus the initial energy density 
can be accessed by a corrected estimation Ec as 



^ = (2A-l)p 



A-l 



dn 



(i?27r)ro dy 



(15) 



Here £b is the Bjorken estimation, which is recovered if 
^ is flat (i.e. A = 1), but for A > 1, both correction fac- 
tors are bigger than 1. Hence the initial energy densities 
are under-estimated by the Bjorken formula. Fig. 3 indi- 
cates fits to BRAHMS pseudo-rapidity distributions from 
ref. [13], these fits indicate that = 8.5 — 10 GeV/fm'^ 
in Au+Au collisions at RHIC. 
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I Normalized pseudorapidity distributions from BRAHMS | 
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FIG. 3: Top panel: Fits of eq. (12) to dn/drj data as mea- 
sured by the BRAHMS collaboration [13] in y^s^ = 200 
GeV Au+Au collisions at various centralities. Bottom panel: 
Sc/bb ratio as a function of r//ro. Using the Bjorken esti- 
mate of £b = 5 GeV/fm^ as given in the BRAHMS White 
Paper [14], and Tf /tq— 7-10, we find an initial energy density 
of £c = (1.7 - 2.0)£s = 8.5 - 10 GeV/fm^ 



IV. LIFE-TIME DETERMINATION 

For a Hwa-Bjorken type of accelerationless, coasting 
longitudinal flow, Sinyukov and Makhlin [15] determined 
the longitudinal length of homogeneity as 



-Rlong — \ —Tb- 

V mt 



Here mt is the transverse mass and tb is the (Bjorken) 
freeze-out time. However, acceleration influences the es- 
timated life-times of the reaction. If the flow is acceler- 
ating, the estimated origin of the trajectories is shifted, 
and we foimd the following correction of eq. (16), at mid- 
rapidity, for broad but finite rapidity distribution: 



Tf Tc 



mt A 



Xtb. 



(17) 



(16) 



Thus in earlier HBT analyses the effective proper-time 
duration of the reaction have been underestimated, as 
pointed out also in refs. [16-19]. BRAHMS pseudo- 
rapidity distributions in Fig. 3 yield A ~ 1.2, and imply 
a 20 % increase in the estimated life-time of the reaction. 



V. SUMMARY 

We have presented new, simple and accelerating solu- 
tions of relativistic hydrodynamics. We have improved 
quantitatively on Bjorken's initial energy density esti- 
mate, by taking into account the longitudinal work, 
which is important for finite rapidity distributions. We 
have fitted the BRAHMS pseudo-rapidity distributions 
with the resulting simple forms, and pointed out that in 
Au+Au coUisions at RHIC, 8.5 - 10 GeV/fm^ initial en- 
ergy densities are reached, a factor of 2 larger, than the 
Bjorken estimate. We have also corrected the Sinyukov- 
Makhlin formula for longitudinal work effects, and found 
an increase of the life-time of the reaction by about 20 %, 
as extracted from the longitudinal HBT radius parameter 
in Au+Au collisions at ^/snn = 200 GeV. 
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